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Abstract 

We consider charged Lepton Flavour Violation (LFV) in the Constrained Minimal 
Supersymmetric Standard Model, extended to include the see-saw mechanism with 
Constrained Sequential Dominance (CSD), where CSD provides a natural see-saw 
explanation of tri-bimaximal neutrino mixing. When charged lepton corrections to 
tri-bimaximal neutrino mixing are included, we discover characteristic correlations 
among the LFV branching ratios, depending on the mass ordering of the right- 
handed neutrinos, with a pronounced dependence on the leptonic mixing angle 
#13 (and in some cases also on the Dirac CP phase 5). 
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1 Introduction 



Over the past decade neutrino physics has revealed the surprising fact not only that 
neutrinos have mass, but also that lepton mixing must involve two large mixing angles, 
commonly referred to as the atmospheric angle 623 and the solar angle #12 [1| . The latest 
neutrino oscillation data [2] is consistent with tri-bimaximal lepton mixing [3]. Theo- 
retical attempts to reproduce this structure typically produce tri-bimaximal mixing in 
the neutrino sector [I], with charged lepton mixing giving important corrections to the 
physical lepton mixing. For example, in the see-saw mechanism [5], sequential domi- 
nance (SD) [6] is well known to provide a natural explanation of hierarchical neutrino 
mass together with large neutrino mixing angles. When certain constraints are imposed 
on the neutrino Yukawa matrix elements then tri-bimaximal neutrino mixing can result 
from such a constrained sequential dominance (CSD) [7j. Charged lepton corrections 
can provide calculable deviations from tri-bimaximal mixing, resulting in predictive neu- 
trino mixing sum rules [3, El [9] which may be proved with future long baseline neutrino 
experiments |10j . 

When neutrino mass models are combined with supersymmetry (SUSY) then lepton 
flavour violation (LFV) is an inevitable consequence [TH [121 in the constrained 
minimal supersymmetric standard model (CMSSM), in which the soft scalar mass ma- 
trices are described by a single universal soft high energy parameter m , and a universal 
trilinear parameter Aq, then the only source of LFV is due to RGE running effects, 
and in this case the connection between LFV processes and neutrino mass models has 
received a good deal of attention [Hj. In the case of SD models it has been shown 
that LFV could reveal direct information about the neutrino Yukawa couplings in the 
diagonal charged lepton basis, depending on the particular nature of the SD, for exam- 
ple whether the dominant right-handed neutrino is the heaviest one or the lightest one 
p~5l [T6] . For example if the dominant right-handed neutrino is the heaviest one, then 
large rates for r — > /ry are expected [151 [IB]- However even in this case, the amount of 
information one can deduce is limited due to the large number of unconstrained Yukawa 
couplings. 

In this paper we consider LFV for the case of CSD, where the number of independent 
neutrino Yukawa couplings is reduced. In this case the LFV predictions are also sensitive 
to the charged lepton mixings, so some further assumptions are required in order to make 
predictions. In addition to tri-bimaximal mixing via CSD, we shall also additionally 
assume CKM-like charged lepton corrections. This will lead to interesting correlations 
in LFV muon and tau decays, independent of the SUSY mass parameters, and Yukawa 
couplings, providing quite specific predictions for LFV. 
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2 (Constrained) Sequential Dominance 



Sequential Dominance (SD) [B] represents classes of neutrino models where large lepton 
mixing angles and small hierarchical neutrino masses can be readily explained within 
the see-saw mechanism. To understand how Sequential Dominance works, we begin by 
writing the right-handed neutrino Majorana mass matrix Mrr in a diagonal basis as 

'M A 

M RR = I M B 1. (1) 



M, 



c , 



We furthermore write the neutrino (Dirac) Yukawa matrix \ v in terms of (1, 3) column 
vectors Ai, B i: Ci as 

Y v = (A B C) , (2) 

using left-right convention. The term for the light neutrino masses in the effective 
Lagrangian (after electroweak symmetry breaking), resulting from integrating out the 
massive right handed neutrinos, is 



-„ ("TA)(Ajv 3 ) {ufB^Bju,) {yFCiKCfvj) 
L '" W A + Mb + W c 



where (i = 1,2,3) are the left-handed neutrino fields. Sequential dominance then 
corresponds to the third term being negligible, the second term subdominant and the 
first term dominant: 

AiAn B«B.\ CjCq , x 

— — > — ^ > . (4) 

M A M B M c y 1 

In addition, we shall shortly see that small #13 and almost maximal #23 require that 

\A X \ < \A 2 \ w |A 3 |. (5) 

Without loss of generality, then, we shall label the dominant right-handed neutrino and 
Yukawa couplings as A, the subdominant ones as B, and the almost decoupled (sub- 
subdominant) ones as C. Note that the mass ordering of right-handed neutrinos is not 
yet specified. Again without loss of generality we shall order the right-handed neutrino 
masses as Mi < M 2 < M 3 , and subsequently identify M A , Mb, Mc with M 1; M 2 , M 3 in 
all possible ways. LFV in some of these classes of SD models has been analysed in 
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Writing A a = \A a \e t(t,Aa , B a = \B a \e l<l>B <*, C Q = \C a \e l<l>c <* and working in the mass 
basis of the charged leptons, under the SD condition Eq. (jl]), we obtain for the neutrino 
mixing angles |6|: 

\M 

tan #23 ~ T~A\ ' ( 6& ) 

I B I 

tan 6\ 2 « J-^ =- , (6b) 

C 2 3\B 2 \ COS fa - S 23 |^3| COS03 

e » „ c ^-^ ) \Bi\(A* 2 B 2 + AlB 3 ) M A e^-^)|A| 

[\A 2 \^ + \A 3 \^ M * + W ( } 

and for the masses 

(W + |4»iy , 

-3 - Wa , (7a) 

I Si IV 



m2 ~ ^HvT' (7b) 

S l2 M B 

mi « 0(|qV/M c ). (7c) 



As in [H] the PMNS phase 5 is fixed by the requirement that we have already imposed 
in Eq. ( 16bl) that tan(#i2) is real and positive, 

c 23 \B 2 1 sin 2 « S23IS3I sin 03 , (8) 

C23IS2I COS02 - S23IS3I COS03 > 0, (9) 

where 

fa = 4>b 2 - <f>B! - <t> + 8 , 

fa = 4>b 3 - <Pb 1 + 4>a 2 - 4>a 3 - <f> + 8 ■ (10) 

The phase is fixed by the requirement (not yet imposed in Eq. fl6cl) ) that the angle #13 
is real and positive. In general this condition is rather complicated since the expression 
for # 13 is a sum of two terms. However if, for example, A\ = then is fixed by: 



-fa-C (11) 

where 

C = &rg (a;b 2 + a;b 3 ). (12) 
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Eq. ffHZ]) may be expressed as 



tan( « 



l£ 2 l 


$23 s 2 + 


-S3 c 23 s 3 


l#2 


S23 c 2 + 


B3 C23C3 



(13) 



Inserting of Eq. (II ip into Eqs. (|H]), (1101) . we obtain a relation which can be expressed 
as 



tan(C + 5) 





B2 C23 s 2 — 


B3 \s 23S 3 




|£ 2 


C23 c 2 + 


-S3 s 23 c 3 



(14) 



In Eqs. ( TT3"j) . (fill) we have written Sj = sin£j, q = cosQ, where we have defined 

C2 = (f>B 2 ~ <pA 2 , C3 = <Pb 3 ~ <pA 3 , (15) 

which are invariant under a charged lepton phase transformation. The reason why the 
see-saw parameters only involve two invariant phases rather than the usual six, is due 
to the SD assumption in Eq. PJ that has the effect of effectively decoupling the right- 
handed neutrino of mass Mc from the see-saw mechanism, which removes three phases, 
together with the further assumption (in this case) of A\ = 0, which removes another 
phase. 



2.1 CSD and tri-bimaximal neutrino mixing 

Tri-bimaximal neutrino mixing [3 J corresponds to the choice [7j: 

(16) 
(17) 

\B S \, (18) 

(19) 

This is called constrained sequential dominance (CSD) [7]. For example, a neutrino 
Yukawa matrix in the notation of Eq. (j2J), which satisfies the CSD conditions in Eqs. ( flBl) - 
( |T9|) . may be taken to be: 



(20) 



1^1 1 


= 0, 


\M 


= \A 3 


\Bx\ 


= \B 2 


A*B 


= 0. 
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where C is not constrained by CSD, since it only gives a sub-sub dominant contribution 
to the neutrino mass matrix, so we have written it as C = (ci, c 2 , c 3 ) above. CSD leads 
to tri-bimaximal mixing in the neutrino mass matrix m„, i.e. to 



fL,tri 



^ijl l/y/3 
-l/VE 1/V3 l/y/2 

l/VE -1/V3 1/V2 



) 



(21) 



3 Charged lepton corrections 



The form of the PMNS matrix will depend on the charged lepton Yukawa matrix whose 
diagonalisation will result in a charged lepton mixing matrix V eL which must be combined 
with Vj} l to form £7pmns- The resulting lepton mixing matrix will therefore not be 
precisely of the tri-bimaximal form, even in theories that predict precise tri-bimaximal 
neutrino mixing. We consider here the case that CSD holds in a basis where the charged 
lepton mass matrix is not exactly diagonal, but corresponds to small mixing. This is a 
situation, often encountered in realistic models [?J [9j [17] . 

In the presence of charged lepton corrections, the prediction of tri-bimaximal neu- 
trino mixing is not directly experimentally accessible. However, this challenge can be 
overcome when we make the additional assumption that the charged lepton mixing 
matrix has a CKM-like structure, in the sense that V eL is dominated by a 1-2 mix- 
ing 6 = 9f 2 , i-e. that its elements (V r eL )i3, (V eh )23, (K L )3i and (V eL )32 are very small 
compared to (V eL )ij = 1,2). In the following, we shall take these elements to be 
approximately zero, i.e. 



where eg = cos 6, sg = sin 6, A is a phase required to diagonalise the charged lepton 
mass matrix [7J, and P is a diagonal matrix of phases P = diag(e mi ,e W2 ,e W3 ) which 
are chosen to remove phases from the product V eL V^ to yield the physical ?7pmns- 111 
the present case it is convenient to choose ui\ — 0, U2 = A, = 0, to yield, 



With this choice, then by constructing L/pmns and comparing to the Standard PDG 
form of this matrix, one obtains, by comparing with Eq.(82) of [TJ, 




(22) 




(23) 



A = 5 - Ti 



(24) 
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where 5 is the Standard PDG CP violating oscillation phase. Also note that A ~ #22 — #12 
where 6ij = argM®-. 

We remark that the assumption that the charged lepton mixing angles are domi- 
nated by (1,2) Cabibbo-like mixing arises in many generic classes of flavour models in 
the context of unified theories of fundamental interactions, where quarks and leptons are 
joined in representations of the unified gauge symmetries [3 [9], [T7] . Under this assump- 
tion, it follows directly from Eq. (1A.3|) that (£/ PM ns)31, (^pmns)32 and (?7 PM ns)33 are 
independent of V eL , and depend only on the diagonalisation matrix VJ^ of the neutrino 
mass matrix. This leads to the parameterization-independent relations [10J: 

|(Vj>i| » |(f/pMNs) 3 i| , (25a) 
\(Vl)s 2 \ « I(t/PMNS) 32 | , (25b) 
1(^)331 « |(E/PMNS)33| ■ (25c) 



In addition to the assumption that V eh is of the form of Eq. (12311 for tri-bimaximal 
neutrino mixing the 1-3 mixing in the neutrino mass matrix is zero, 

(KJ L )i3 = 0. (26) 

Using Eq. ( 1261) and applying the standard PDG parameterization of the PMNS matrix 
(see e.g. p]), Eq. (l2Hal) leads to the sum rule El E]: 

S 23 S 12 ~ | s 23-Sl2 - Sl3C 2 3Cl 2 e l5 | « S 23 Sl2 - S 13 C 23 C 12 COs(5) , (27) 

where the last step holds to leading order in S13. This sum rule can be used to test 
tri-bimaximal [d\ 2 = arcsin(-^g)) structure of the neutrino mass matrix in the presence 
of CKM-like charged lepton corrections. 



4 LFV in CSD with charged lepton corrections 

When dealing with LFV it is convenient to work in the basis where the charged lep- 
ton mass matrix is diagonal. Let us now discuss the consequences of charged lepton 
corrections of the form of Eq. (T23]) for the neutrino Yukawa matrix with CSD. After re- 
diagonalising the charged lepton mass matrix, resulting in the assumed charged lepton 
mixing matrix in Eq. (I2"3l . Y v in Eq. (|2"U]) becomes transformed as: 

Y v - Y' v = V eL Y v . (28) 
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In the diagonal charged lepton mass basis the neutrino Yukawa matrix therefore be- 
comes: 

Sge -i\ e ifo & fa _ See -^) e ifr ( ClQ? _ c 2 s e e- iX ) 
Yl=(A' B' C')=\ -acge^ 3 b {c e + s e iX )e^ ( Cl s e e iX + c 2 c e ) |. (29) 

ae % be if)2 c 3 

where the column vectors A', B', C are now defined in the diagonal charged lepton 
basis according to Eq. fl2"9"]) . Thus the results in Eqs. © with the redefined column 
vectors A', B', C now yield the physical lepton mixing angles since these are equal to 
the neutrino mixing angles in the diagonal charged lepton basis of Eq. ( 1291 ) . 

At leading order in a mass insertion (MI) approximation [TT1 [T2] the branching 
fractions of LFV processes are given by 

a 3 

B Tij = Br(k -> Iff) « -^f(M 2 , //, m e ) |m|.. tan 2 (3 , (30) 

where l\ — e, l 2 — £3 = r, and where the off-diagonal slepton doublet mass squared is 
given in the leading log approximation (LLA) of the CMSSM by 

mf "> « -<M±^)^ , (31) 

with the leading log coefficients given by 

K 21 = A ^A'*\n^- + B' 2 B'*\n^- + C' 2 C'*\n A 



M A z 1 "' M B ' ~' Z ~ L "' M c ' 

A' 31 = ^rinA + ^ Bri „^ + CM -lnA. (32) 

The factors ^ in Eq. ( 1301) represent the ratio of the leptonic partial width to the total 
width, 

_ r(z f -> fji/^j) 
™~ r(/^aii) • 

Clearly £21 = 1 but £32 is non-zero and must be included for correct comparison with 
the experimental limit on the branching ratio for r — ► ^7. This factor is frequently 
forgotten in the theoretical literature. 
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If LFV is only induced by RG effects from Y' v on the soft breaking terms, as in the 
CMSSM, then in the LLog and MI approximation, the branching ratios for LFV charged 
lepton decays, like l{ — > £jj, are proportional to 

Br l3 oc \K l3 \ 2 = \(A'A\ J HA/M A ) + (B'B'^ l ,\n(A/M B ) + (CV\ 3 \n(A/M c )\ 2 \34) 

We have only assumed so far that the right-handed neutrino mass matrix has the diago- 
nal form shown in Eq. (OQ), Mrr = diag(M4, Mg, Mc) with the dominant right-handed 
neutrino labelled by A, the leading subdominant one labelled by B, and the decoupled 
one labelled by C. However the masses of the right-handed neutrinos are not yet or- 
dered, and we have not yet specified which one is the lightest and so on. After ordering 
Ma, Mb, Mc according to their size, there are six possible forms of Y£ obtained from 
permuting the columns, with the convention always being that the dominant one is la- 
belled by A, and so on. In particular the third column of the neutrino Yukawa matrix 
could be A', B' or C depending on which of Ma, M b or Mc is the heaviest. 

In hierarchical models, the (3, 3) elements of the Yukawa matrices describing quarks 
and charged leptons are amongst the largest elements in the Yukawa matrices. In unified 
models this will also be the case for the neutrino Yukawa matrix. If the heaviest right- 
handed neutrino mass is Ma then the third column of the neutrino Yukawa matrix will 
consist of the A 1 column, and since Y 33 = A' 3 and A\ ~ A' 2 ~ A' 3 ~ a then we conclude 
that all elements of A' must dominate over those of B', C, and hence all LFV processes 
will be determined approximately by (A'A*)ij. Similarly if the heaviest right-handed 
neutrino mass is Mb then the third column of the neutrino Yukawa matrix will consist 
of the B' column, and since Y 33 = B 3 and B[ ~ B' 2 ~ B 3 ~ b then we conclude that 
all elements of B' must dominate over those of A', C, and hence all LFV processes 
will be determined approximately by (B'B^)ij. Finally if the heaviest right-handed 
neutrino mass is Mc then the third column of the neutrino Yukawa matrix will consist 
of the C column which contains the large element Y 33 = C 3 . However in this case 
we cannot conclude that all elements of C must dominate over those of A', B' for the 
determination of LFV processes since the elements c\, C2 are undetermined by the see- 
saw mechanism and could even be set equal to zero. Nevertheless it is possible that 
in this case all elements of C could dominate over those of A', B' and hence all LFV 
processes could be determined approximately by (C'C t)^. In the following we consider 
the LFV predictions arising from the three cases 

M 3 = M A , M 3 = Mb , M 3 = M C , (35) 

corresponding to the dominant Yukawa columns being A', B', C, respectively. 
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5 Predictions for the ratios of LFV branching ratios 



After ordering Ma, Mb, Mq according to their size, there are six possible forms of Y' v 
obtained from permuting the columns, with the convention always being that the dom- 
inant one is labeled by A', and so on. In particular the third column of the neutrino 
Yukawa matrix could be A', B' or C depending on which of Ma, Mb or Mq is the 
heaviest. If the heaviest right-handed neutrino mass is Ma then the third column of 
the neutrino Yukawa matrix will consist of the (re-ordered) first column of Eq. (1291) and 
assuming Y 3 % ~ 1 we conclude that all LFV processes will be determined approximately 
by the first column of Eq. ( 1291) . Similarly if the heaviest right-handed neutrino mass is 
Mb then we conclude that all LFV processes will be determined approximately by the 
second column of Eq. (129"]) . Note that in both cases the ratios of branching ratios are 
independent of the unknown Yukawa couplings which cancel, and only depend on the 
charged lepton angle 9 = 9f 2 (and in some cases on A) , which in the case of tri-bimaximal 
neutrino mixing is related to the physical reactor angle by 6*13 = 0f 2 / v2 = 9 / ' s/2 (7J [9]. 
Also note that A = 5 — 71 where 5 is the Standard PDG CP violating oscillation phase. 
The predictions for these two cases will now be discussed in detail. We will also com- 
ment on the third case M 3 = Mc, which is less predictive, and give an explicit minimal 
example. 



5.1 M 3 = M A 

In this case, assuming that the third column of the neutrino Yukawa matrix (associated 
with the heaviest right-handed neutrino and hence the largest Yukawa couplings) is the 
dominant column A' associated with the atmospheric neutrino of mass 777.3, one can read 
off from Eq. flM} and Eq. fl2HD that the Br^ = Br(£i -> ^7) now satisfy 

Br^ oc \a 2 s e c e \ 2 , (36) 
Br re oc \a 2 s e \ 2 ire, (37) 
Br Tfl oc \a 2 \ 2 ^. (38) 

Note that 9 = 0\ 2 = V%9i3, so there is a direct (and simple) connection to the measurable 
lepton mixing angle #13 in neutrino oscillation experiments in this case. In particular, 
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Figure 1: The left panel shows the ratios of branching ratios Br^ of LFV processes ii — > £j7 in CSD 
for M 3 = M A with right-handed neutrino masses Mi = 10 8 GeV, M 2 = 5 x 10 8 GeV and M 3 = 10 14 
GeV. Here the solid lines show the (naive) prediction, from the MI and LLog approximation and with 
RG running effects for the other parameters neglected, while the dots show the explicit numerical 
computation (using SPhcno2.2.2. [T!5] extended by software packages for LFV branching ratios and 
neutrino mass matrix running [501 HI]) with universal CMSSM parameters chosen as mo = 750 GeV, 
mi/2 — 750 GeV, A$ = GeV, tan/? = 10 and sign(/i) = +1. The right panel shows the predictions 
(from full computation) for Br^ e = Br(n — > &y) in the CMSSM extended by the see-saw mechanism 
with CSD for the case M 3 = M A with 13 = 3° and 5 = 0. In this panel we have chosen the CMSSM 
parameters to satisfy Aq — GeV, tan/3 = 10 and TO1/2 = 5mo, which approximately corresponds to 
the successful stau co-annihilation region of LSP neutralino dark matter (DM) giving CIdm within the 
current WMAP limits. 



we predict 
Br 



fie 



Br 



Br 



Br Tl 
Br Tt 
Br T . 



1 2 S^te 



(*) 



2^e 



I sin(2v^0i 3 
2 t 



cos(\/2^ 
sin(V20 13 ) 



Lie 



(39) 
(40) 
(41) 



The predictions for the ratios of branching ratios as a function of #13 as well as for Br^ e , 
for some sample choice of parameters, are shown in Fig. [H 



10 



5.2 M 3 = M B 



In this case, assuming that the third column of the neutrino Yukawa matrix (associated 
with the heaviest right-handed neutrino and hence the largest Yukawa couplings) is the 
leading subdominant column B' associated with the solar neutrino of mass m2, one can 
read off from Eq. ( 1341) and Eq. (1291) that the Br^ = Br(£i — > ^7) now satisfy 



Br 



Br Te 
Br 



oc \b 2 (c e - s e e~ iX )(c e + s e e iX ) 
oc \b 2 (c e 



lie 1 



s s e^ x )\ 2 U, 



oc \b 2 (c e + s e e lX )\ 2 ^. 



(42) 
(43) 
(44) 



bmce 9 = 9 e 12 = V29 13 , there is again a connection to the measurable lepton mixing 
angle #13 in neutrino oscillation experiments. Furthermore, the branching ratios also 
depend on the phase A, which is related to the Standard PDG CP violating oscillation 
phase 6 by A = 5 — 71. The ratios of branching ratios are predicted as 



Br TIJj 


c g -s e e iX \ 


2 S/xe 


-- | cos(v / 2^i 3 ) + sin(^20 13 )e- 


-iS 2 S/xe 

STjtl 


Br^ 
Br Te 




| cos(V20 13 ) - sm(V26 n )e iS 


2^e 


Br Te 


c e - s e e iX 
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cos(V26 13 ) + sm(V26 13 )e- iS 
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Br Tfl 


ce + s e e iX 




cos(V2^i 3 ) - sin(V2^i 3 )e i<5 





(45) 
(46) 

(47) 



Fig. [2] shows the predictions for the ratios of branching ratios as a function of #13, for the 
example 5 = 0, as well as the prediction for Br^ e for some sample choice of parameters. 



5.3 M 3 = M c 

In this case, assuming that the third column of the neutrino Yukawa matrix (associated 
with the heaviest right-handed neutrino and hence the largest Yukawa couplings) is the 
most subdominant column C associated with the lightest neutrino of mass mj, assuming 
that c 3 ~ 1, one can see from Eq. flMl) and Eq. fl2TJ|) that the Br^ = Br(£i — * ^7) now 
depend on undetermined coefficients C\,C2- Hence we cannot make definite predictions. 
Moreover, in some cases, the subdominant column of Yukawa coupling also contributes 
at the same order as the dominant one. Nevertheless, charged lepton corrections also 
have an impact here. Let us therefore generalize V eL to include also a small 8 23 <C 9f 2 . 
As a minimal case, let us furthermore consider 

C = (0,0, c) T . (48) 
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013 [°] log 10 (M 3 /GeV) 

Figure 2: The left panel shows the ratios of branching ratios Brij of LFV processes £i — > ^7 in CSD 
for M3 = Mb and for 5 = 0. The other parameters are chosen as in Fig. [T] The solid lines in the 
left panel show the (naive) prediction, from the MI and LLog approximation and with RG running 
effects neglected, while the dots show the explicit numerical computation. The right panel shows the 
predictions (from full computation) for Br^ e — Br(/i — > ej) in the CMSSM extended by the see-saw 
mechanism with CSD for the case M3 = Mb with #13 = 3° and <5 = 0. The other parameters are chosen 
as in Fig. [TJ 



This may be viewed as minimal scenario regarding LFV, since typically (barring can- 
cellations) the zeros are replaced by small entries and since, as mentioned above, the 
subdominant column of Y v can not in general be neglected. For a more accurate treat- 
ment of this scenario with respect to the charged lepton corrections, the (typically) even 
smaller 9f 3 -C 9 23 can be included analogously. 

Including charged lepton corrections from 9\ 2 and 9 23 (by Y u — > V ei Y v ) leads to 
approximately 

C = (cs e 23 s e 12} cs e 23 ,c) T , (49) 
and thus to the following relations for the branching ratios: 

Br^ oc |c 2 (^ 3 ) 2 ^ 2 | 2 ^ e , (50) 

Br re OC |c S23 S 12l £i~e > (51) 

Br Ti , cc \c 2 s e 23 \ 2 i Ttl . (52) 

As in the cases M 3 = Ma and M 3 = M^, the relation 6\ 2 = \/29i 3 holds under the con- 
sidered assumption about the charged lepton corrections. For the ratios of the branching 
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2 4 6 8 10 12.5 13 13.5 14 14.5 15 15.5 
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Figure 3: The left panel shows the ratios of branching ratios Br.y of LFV processes li — > £jj for a 
minimal example with CSD and M3 = Mc described in the text. The solid lines in the left panel show 
the (naive) prediction, from the MI and LLog approximation and with RG running effects neglected, 
while the dots show the explicit numerical computation. The right panel shows the predictions (from 
full computation) for Br^ e — Br(fi — > ej) in the CMSSM extended by the see-saw mechanism with 
CSD for the case M3 = Mc in the scenario with mi = 10 -3 eV and 6 = 0. The other parameters are 
chosen as in Fig. Q] 



ratios we obtain 

2 ^ , (53) 

(54) 
(55) 

The predictions for the ratios of branching ratios as a function of #13 as well as for Br^ 
as a function of #13 and mm (set equal to 5m as an example) are shown in Fig. [3] To 
give an explicit example, we have chosen s| 3 = sin(#^ 3 KM ) « 2.36° and other parameters 
as stated in the caption of Fig. [TJ We would like to stress again that, in contrast to 
the cases M 3 = Ma and M 3 = Mb discussed above, the shown results are no definite 
predictions for the case M3 = Mc, but rather order of magnitude examples for certain 
classes of models of CSD where the LFV branching ratios are significantly smaller than 
for CSD with M3 = Ma and M3 = Mb- As can be seen from Fig. [3j this scenario can 
be readily distinguished from the cases M 3 = Ma and M 3 = M B . 
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6 Conclusions 



We have considered charged Lepton Flavour Violation (LFV) in the Constrained Mini- 
mal Supersymmetric Standard Model, extended to include the see-saw mechanism with 
Constrained Sequential Dominance (CSD), where CSD provides a natural see-saw expla- 
nation of tri-bimaximal neutrino mixing. When Cabibbo-like charged lepton corrections 
to tri-bimaximal neutrino mixing are included, this leads to characteristic correlations 
among the LFV branching ratios Br Tfl , Br^ e and Br Te which may be tested in future 
experiments. 

There are two main differences between the study here and that in [16] where pre- 
dictions for LFV were also presented for the CMSSM with SD. The first difference is 
that here we have focused on the special case of CSD, corresponding to tri-bimaximal 
neutrino mixing, where the neutrino Yukawa couplings are very tightly constrained com- 
pared to the general SD case. The second difference is that we have considered the effect 
of charged lepton corrections, which were not included in [IB] . In particular we have 
mainly considered Cabibbo-like charged lepton corrections, which when combined with 
CSD leads to a very tightly constrained scenario in which ratios of branching ratios 
depend on # 13 , which is related to the charged lepton mixing angle 9f 2 - The predic- 
tions also depend crucially on which column of the Yukawa matrix is associated with 
the heaviest right-handed neutrino M 3 , since this column will have the largest Yukawa 
couplings. 

For the case M3 = Ma, also known as Heavy Sequential Dominance (HSD) since 
the dominant right-handed neutrino is the heaviest one, we find the characteristic ratios 
in Fig. 1. Compared to the results in [16], the hierarchy between Br^ e and Br T/1 is 
much milder. This can be understood from the fact that in [IB] it was assumed that 
I Ai| <C |A 2 | « I A 3 1 ~ 1 (ignoring charged lepton corrections) which led to large Br Tfl but 
small Br^ e . However, including charged lepton corrections, we see that \A X \ ~ \A' 2 \ ~ 
I A 3 1 ~ 1, leading to both large Br Til and large Br^ e . In the present case we focus on tri- 
bimaximal neutrino mixing, which before charged lepton corrections are included implies 
that |Ai| <C |A 2 | = I A3 1 ~ 1, corresponding to the CSD explanation of tri-bimaximal 
neutrino mixing. Then, after Cabibbo-like charged lepton corrections are included, this 
leads to well defined predictions for the each of the couplings \A[\, |A' 2 |, |A' 3 |, and hence 
rather precise predictions for ratios of Br Tfl , Br^ e and Br Te , which depend on #13, as 
shown in Fig. 1. We reemphasize that, after charged lepton corrections are included, 
I Ai| ~ |A' 2 | ~ I A 3 1 ~ 1, and hence both Br Tfl and Br^ are large in this case, unlike [16] 
where charged lepton corrections were ignored. 

In the case M 3 = Mb, where the leading subdominant right-handed neutrino respon- 
sible for the solar neutrino mass is the heaviest one, the predicted ratios of branching 
ratios are even milder, corresponding to the fact that all the Yukawa coupling in this col- 
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umn are equal before the inclusion of charged lepton corrections, \Bi \ = \B 2 \ = \B 3 \ ~ 1, 
again corresponding to the CSD explanation of tri-bimaximal neutrino mixing. When 
Cabibbo-like charged lepton corrections are included this again leads to characteristic 
predictions for ratios of Br T/I , Br^ e and Br Te , depending on #13 (shown in Fig. 2 for 
5 = 0) as well as on the Dirac CP phase 5, as given in Eqs. (j4"5l) - ( 1471) . 

The least predictive case is M3 = Mc, which includes the case where the dominant 
right-handed neutrino is the lightest one known as light sequential dominance (LSD). In 
this case the generic prediction from [TB] was that the Br Tll was generally quite small, 
typically of order Br^ e , due to the small neutrino Yukawa couplings. In particular the 
neutrino Yukawa couplings of the third column were not considered relevant due to the 
large mass of the associated right-handed neutrino, which was assumed to exceed the 
GUT scale from which the RGEs were run down. Then the relevant Yukawa couplings 
were those from the second column, which all take similar (small) values leading to 
Br TfJL ~ Br^. This may also be the case here, since including charged lepton corrections 
will not change this result, and CSD will only strengthen this conclusion. However in 
other cases, for example if the RGEs are run from the Planck scale, the third column of 
the neutrino Yukawa matrix should not be ignored. In Fig. 3 we considered an example 
of this, in which the LFV arises solely from the third column, and the Yukawa couplings 
in this column are again determined from charged lepton corrections, assuming that 
C = (0, 0, c) T , which may be approximately true in practice, but which is by no means 
guaranteed. 

In summary, the results presented here once again confirm that Br Tfl , Br^ e and 
Br Te are all expected to be observed in the (near) future. If they are observed in the 
ratios predicted here, for some value of #13, then this may be an indication of a high 
energy theory with the characteristics of the CMSSM extended to include the see-saw 
mechanism with CSD, corresponding to tri-bimaximal neutrino mixing corrected by 
Cabibbo-like charged lepton mixing angles. 

Acknowledgements 

S. A. would like to thank M. J. Herrero and E. Arganda for many very helpful discussions 
about LFV in SUSY see-saw scenarios and for providing their software packages on LFV 
muon and tau decays. 



15 



Appendix 



A Conventions 

In general, the mixing matrix in the lepton sector, the PMNS matrix £7pmnS; is defined 
as the matrix which appears in the electroweak coupling to the W bosons expressed in 
terms of lepton mass eigenstates. With the mass matrices of charged leptons M e and 
neutrinos m u written as 

£ = -e L M e e R - \v L m v vl + H.c. , (A.l) 
and performing the transformation from flavour to mass basis by 

V eh M e Vj R = diag(m e , m M , m T ), V UL m u Vj L = diag(mi, m 2 , m 3 ), (A. 2) 

the PMNS matrix is given by 

f/pMNS = V eL Vl. (A.3) 

Here it is assumed implicitly that unphysical phases are removed by field redefinitions, 
and t/pMNS contains one Dirac phase and two Majorana phases. The latter are physical 
only in the case of Majorana neutrinos, for Dirac neutrinos the two Majorana phases 
can be absorbed as well. 

The standard PDG parameterization of the PMNS matrix (see e.g. [18J ) is: 

C12C13 S12C13 s 13 e~ lS 

UpMNS = I -C23S12 - Si 3 S23Ci2e 1<5 C23C12 - S 13 S 2 3Sl 2 e l5 S23C13 I -Pjvlaj , (A.4) 
S23S12 - Si 3 C23Ci 2 e 1<5 -S23C12 - Sl3C 2 3Sl2e l5 C23C13 

which is used in most analyses of neutrino oscillation experiments. Here S is the so- 
called Dirac CP violating phase which is in principle measurable in neutrino oscillation 

. OCX .02 

experiments, and PMaj = diag(e' 2 , e 1 2 , 0) contains the Majorana phases a±, a%. In the 
following we will use this standard parameterization (including additional phases) also 
for Vj L and denote the corresponding mixing angles by 6^, while the mixing angles 
without superscript refer to the PMNS matrix. 
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